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Abstract. In this note we present the main results from the recent work [5], 
which for the first time establish Landau damping in a nonlinear context. 

Keywords. Landau damping; plasma physics; astrophysics; Vlasov-Poisson 
equation. 

1. Introduction 

The "standard model" of classical plasma physics is the Vlasov-Poisson-Landau 
equation [T3], E] , here written with periodic boundary conditions and in adimensional 
units: 

d f log A 

(1) -£ + v- VJ + F[f] ■ VJ = Q L (f, f), 

where / = f(t,x,v) is the electron distribution function (t > 0, v G M 3 , x G T 3 = 

R 3 /Z 3 ), 

(2) F\f]{t,x) = - jj VW(x-y)f(t,y,w)dwdy 

is the self-induced force, W(x) = l/\x\ is the Coulomb interaction potential, and Ql 
is the Landau collision operator, described for instance in [12]. The parameter A is 
very large, ranging typically from 10 2 to 10 30 . 

On very large time scales (say O (A/ log A)), dissipative phenomena play a non- 
negligeable role, and the entropy increase is supposed to force the (slow) convergence 
to a Maxwellian distribution. Thanks to the recent progress on hypocoercivity, this 
mechanism is now rather well understood, as soon as global smoothness estimates 
are available (see [UJ and the references therein). 

Ten years after devising this collisional scenario, Landau [B] formulated a much 
more subtle prediction: the stability of homogeneous equilibria satisfying certain 
conditions — for instance any function of \v\, not necessarily Gaussian — on much 
shorter time scales (say 0(1)), by means of purely conservative mechanisms. This 
phenomenon, called Landau damping, is a property of the (collisionless) Vlasov 
equation, obtained by setting A = oo in ([T]). This is a theoretical cornerstone of 
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the classical plasma physics (among a large number of references let us mention [T]). 
Similar damping phenomena also occur in other domains of physics. 

The Landau damping has been since long understood at the linearized level [3], 
[10], but the study of the full (nonlinear) equation poses important conceptual 
and technical problems. As a consequence, up to now the only existing results were 
proving existence of some damped solutions with prescribed behavior as t — > ±00 
[21 [5]. We fill this gap in a recent work [9], whose main result we shall now describe. 

2. Main result 

If / is a function defined on T d x M. d , we note, for any fc£Z d and 77 G Mr, 
f(k,v) = I f(x,v)e- 2mk - x dx, f{k,rj)= If f(x, v) e" 2 ^ e~ 2l ^ v dv dx. 



(3) ll/IU,^ = sup(|/(M)| ^ XH ^) +11 I/O*. «)l ^ H dvdx 



We also set, for A, /x, /5 > 0, 

'(' , ,, 

Theorem 1 (nonlinear Landau damping for general interaction). Let d > 1, and 
f° : M. d — > M + an analytic velocity profile. Let W : T d — > R &e an interaction 
potential. For any k G Z d ; ^ G C, we set 

/"OO 

£(fc,f) = -4vr 2 l?(A;) / e 2 ^**|/°(A;t)||A;| 2 tdt. 

JO 

We assume that there is A > suc/i i/ia£, for e small enough, 

(4) sup |/°(t?)| e 2 ^l < Co, V -rllV^l^^) < L7 , 

(5) inf inf |£(fc,0 - 1| > k > 

kez d 0<3?C<A 1 1 



(6) 3 7 >l;VA;GZ d ; |W(Ar)| < 

Then as soon as < A' < A, < // < /3 > 0, r G N, taere are £ > and C > 0, 
depending on d, 7, A, A', /1, //, Co, ft, Ciy, f3, r, such that if fa > satisfies 

(7) 5:= ||/ 4 -/°|U^<e, 
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then the unique solution of the nonlinear Vlasov equation 

(8) ^+v-V x f+F[f]-V v f = 0, F[f](t,x) = - J J VW(x-y)f(t,y,w)dwdy, 

defined for all times and such that /(O, ■) = f%, satisfies 

(9) \\p(t, O-Pooll^^C^e- 2 ^, 

where p(t,x) = J f(t,x,v)dv,p ao = J J fi(x,v) dv dx. Futhermore, there are analytic 
profiles f +OQ (v), /_oo(t>) such that 

f(t, ■ ) ^ ±0 °> / ±00 weakly 

f{t, x, ■ ) dx W± °°> f ±ao strongly (in C r (R d v )), 

these convergences being also 0(5 e~ 27rA ''*'). 

This theorem, entirely constructive, is almost optimal, as the following comments 
show. 

Comments on the assumptions: The periodic boundary conditions of course are 
debatable; in any case, the counterexamples of Glassey and Schaeffer [I] show that 
some confinement mechanism — or at least a limitation on the wavelength — is 
mandatory. Condition (jlj) quantitatively expresses the analyticity of the profile f°, 
without which we could not hope for an exponential convergence. The inequality 
(jSJ) is a linear stability condition, roughly optimal, covering all physically interesting 
cases: in particular the (attractive) Newton interaction for wavelengths shorter than 
the Jeans unstability length; and the (repulsive) Coulomb interaction around radially 
symmetric profiles /°, for all wavelengths. On the other hand, condition ([6]) shows 
up only in the nonlinear stability; it is satisfied by Coulomb and Newton interactions 
as a limit case. As for the condition (j7|), its perturbative nature is natural in view 
of theoretical speculations and numerical studies in the subject. 

Comments on the conclusions: 

(1) The large-time convergence is based on a reversible, purely deterministic 
mechanism, without any Lyapunov functional neither variational interpreta- 
tion. The asymptotic profiles f±oo eventually keep the memory of the initial 
datum and the interaction. This convergence "for no reason" was not really 
expected, since the quasilinear theory of Landau damping [U Vol. II, Section 
9.1.2] predicts convergence only after taking average on statistical ensembles. 



4 



C. MOUHOT AND C. VILLANI 



(2) This result can be interpreted in the spirit of the KAM theorem: for the 
linear Vlasov equation, convergence is forced by an infinity of invariant sub- 
spaces, which make the model "completely integrable" ; as soon as one adds 
a nonlinear coupling, the invariance goes away but the convergence remains. 

(3) Given a stable equilibrium profile f°, we see that an entire neighborhood - 
in analytic topology — of f° is filled by homoclinic or (in general) heteroclinic 
trajectories. Only infinite dimension allows this remarkable behavior of the 
nonlinear Vlasov equation. 

(4) The large time convergence of the distribution function holds only in the 
weak sense; the norms of velocity derivatives grow quickly in large time, 
which reflects a filamentation in phase space, and a transfer of energy (or 
information) from low to high frequencies ("weak turbulence"). 

(5) It is this transfer of information to small scales which allows to reconcile the 
reversibility of the Vlasov-Poisson equation with the seemingly irreversible 
large-time behavior. Let us note that the "dual" mechanism of transfer of 
energy to large scales, also called radiation, was extensively studied in the 
setting of Hamiltonian systems. 

Much more comments, both from the mathematical and the physical sides, can 
be found in [9]. 



3. Linear stability 

The linear stability is the first step of our study; it only requires a reduced technical 
investment. 

After linearization around a homogeneous equilibrium f°, the Vlasov equation 
becomes 

(10) — + v-V x h- (W * p) ■ V,/° = 0, p = J halv. 

It is well-known that this equation decouples into an infinite number of independent 
equations governing the modes of p: for all k G Z d and t > 0, 

(11) p(t,k)- f K°(t-T,k)p(r,k)dT = h,(k,kt), 

Jo 

where hi is the initial datum, and K° an integral kernel depending on f°: 

(12) K°(t,k) = -A7r 2 W{k)f°{kt) \k\ 2 t. 
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Then from classical results on Volterra equations we deduce that for all k 7^ the 
decay of p(t, k) as t — > 00 is essentially controlled by the worst of two convergence 
rates: 

• the convergence rate of the source term in the right-hand side of pip , which 
depends only on the regularity of the initial datum in the velocity variable; 

• e~ xt , where A is the largest positive real number such that the Fourier-Laplace 
transform (in the t variable) of K° does not approach the value 1 in the strip 
{0 < ^tz < A} C C. The problem lies in finding sufficient conditions on f° to 
guarantee the strict positivity of A. 

Since Landau, this study is traditionally performed thanks to the Laplace trans- 
form inversion formula; however, with a view to the nonlinear study, we prefer a 
more elementary and constructive approach, based on the plain Fourier inversion 
formula. 

With this method we establish the linear Landau damping, under conditions (jSJ) 
and (jlj), for any interaction W such that \7W £ L 1 (T d ), and any analytical initial 
condition (without any size restriction in this linear context). We recover as par- 
ticular cases all the results previously established on the linear Landau damping 
[31 El EE] ; but we also cover for instance Newton interaction. Indeed, condition (J3J) 
is satisfied as soon as any one of the following conditions is satisfied: 

(a) VA; £ Z d , Vz £ M, W{k) > 0, z<j/ k (z) < 0, where <p k is the "marginal" of f° 
along the direction k, defined by 



To establish the nonlinear stability, we start by introducing analytic norms which 
are "hybrid" (based on the size of derivatives in the velocity variable, and on the size 
of Fourier coefficients in the position variable) and "gliding" (the norm will change 
with time to take into account the transfer to small velocity scales). Five indices 





We refer to [HI Section 3] for more details. 
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provide all the necessary flexibility: 
(13) 



LP(dv) 



J2 (1 + W (V„ + 2t7rrk) n f(k 1 

keZ d nGN d 

(By default 7 = 0.) A tedious injection theorem "a la Sobolev" compares these 
norms to more traditional ones, such as the H/Ha,^,/? norms appearing in (j3J). 

The Z norms enjoy remarkable properties with respect to composition and prod- 
uct. The parameter r partly compensates for filament at ion. Finally, the hybrid 
nature of these norms is well adapted to the geometry of the problem. If / depends 
only on x, the norm (fT3|) coincides with the norm jF Ar+At ' 7 defined by 

(14) ||/||^*, = E e 2w(A ^ )|fc| (1 + \k\y. 

(We also use the "homogeneous" version jF Ar+/i ' 7 where the mode k = is removed.) 

Then the Vlasov equation is solved by a Newton scheme, whose first step is the 
solution of the linearized equation around /°: 

r = f + h 1 + . . . + h n , 



(15) 

(16) 
n > 1, 




'd t h l + v ■ S7 x h l + Flh 1 } ■ V v f° = 
h\0,-) = f l -f 

v ■ V x h n+l + F[f n ] ■ V v h n+1 + F[h n+l ] ■ 
0. 



WJ n = -F[h n ] ■ V v h r 



In a first step, we establish the short-time analytic regularity of /i n (r, •) in the 
norm Zr X \ this step, in the spirit of a Cauchy-Kowalevskaya theorem, is per- 
formed thanks to the identity 



;i7) 



dt 



ll/l 



K 



1+T 



liv/ll 



^A(t),m(t);p , 



t=T 

where X(t) = A — Kt, fi(t) — fi — Kt. 

In a second step, we establish uniform in time estimates on each h n , now with a 
partly Eulerian and partly Lagrangian method, integrating the equation along the 
characteristics {X™ t , V™ t ) created by the force F[f n ]. (Here r is the initial time, 
t the current time, (x, v ) the initial conditions, (X n ,V n ) the current conditions.) 
The smoothness of these characteristics is expressed by controls in hybrid norm on 
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the operators Q™ T (x, v) = (X" T , V t n T )(x + v(t — r),v), which compare the perturbed 
dynamics to the unperturbed one; these are informally called scattering operators. 

Then we propagate a number of estimates along the scheme; the most important 
are (slightly simplifying) 



sup 

T>0 



h n (r, • , v) dv 



<S r , 



(19) 



(20) 



sup 

t>T>0 



B 



7 A n (l + b),/j 



nil < 5 n , b = bit) = 



? A n (l + b),(/j n ,7);oo 

't—K- 
1+6 



^ E 



4e' 



-27r(A fe -A n+ i)r 



,fc=l 



27r(A fc - X n+i y 



min{i — r ; 1}. 



Notice, in ffTS]) . the linear increase in the regularity of the spatial density, which 
comes at the same time as the deterioration of regularity in the v variable. In (fT9~l) . 
the additional time-shift in the indices by the function b(t) will be crucial to absorb 
error terms coming from the composition; the constant B itself is determined by the 
previous small-time estimates. Finally, in (120| . notice the uniform in t control, and 
the improved estimates in the limit cases t — > r and r — > oo; also this is important 
for handling error terms. The constants A n and \x n decrease at each stage of the 
scheme, converging — not too fast — to positive limits Aoo, /i^; at the same time, 
the constants S n converge extremely fast to 0, which guarantees "by retroaction" the 
uniformity of the constants in the right-hand side of (1201) . 

The estimates (l20l are obtained by repeated application of fixed point theorems 
in analytic norms. Another crucial ingredient to go from stage n to stage n + 1 is the 
mechanism of regularity extortion, which we shall now describe in a simplified 
version. Given two distribution functions / and /, depending on t, x, v, let us define 



a(t, x) 



[ [ {F[f}-V v f){r,x-v(t-r),v)dvdT. 
Jo J«. d 



This quantity can be interpreted as follows: if particles distributed according to 
/ exert a force on particles distributed according to /, then a is the variation of 
density J f dv caused by the reaction of / on /. We show that if / has a high 
gliding regularity, then the regularity of a in large time is better than what would 
be expected: 



(21) 



Ht, -)IU+m< [ K(t,r)\\F[f(, 
Jo 



I jpAr + /j,7 



dr. 
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where 

'|VJ( S , 



K(t,r) 



sup 

0<s<t 



e -2ir(\-\)\k(t-T)+£r\ p -27r(/Z-/i)|^| 

1+r) sup 



life 



The kernel K(t,r) has integral 0(t) as t — > oo, which would let us fear a violent 
unstability; but it is also more and more concentrated on discrete times r = kt/(k — 
£); this is the effect of plasma echoes, discovered and experimentally observed in 
the sixties [7]. The stabilizing role of the echo phenomenon, related to the Landau 
damping, is uncovered in our study. 

Then we analyze the nonlinear response due to echoes. If 7 > 1, from (|2~Tj) one 
deduces that the response is subexponential, and therefore can be controlled by an 
arbitrarily small loss of gliding regularity, at the price of a gigantic constant, which 
later will be absorbed by the ultrafast convergence of the Newton scheme. In the 
end, part of the gliding regularity of / has been converted into a large-time decay. 

When 7 = 1, a finer strategy is needed. To handle this case, we work on the 
response mode by mode, that is, estimating the size of f)(t, k) for all k, via an 
infinite system of inequalities. Then we are able to take advantage of the fact that 
echoes occurring at different frequencies are asymptotically rather well separated. 
For instance, in dimension 1, the dominant echo occurring at time t and frequency 
k corresponds to r = kt/(k + 1). 

In practice, straight trajectories in (T2"Tj) must be replaced by characteristics (this 
reflects the fact that / also exerts a force on /), which is a source of considerable 
technical difficulties. Among the tools used to overcome them, let us mention a 
second mechanism of regularity extortion, acting in short time and close in spirit to 
velocity- averaging lemmas; here is a simplified version of it: 

(22) \\a(t, OII^+m < / \\ F [f( T ' • )] HjTAiT-Kt-rji+ivy ||V/(t, ^W^a+t),^):! dr. 

JO T-bt/(l+b) 

We see in (|2"2"|) that the regularity of a is better than that of F[f], with a gain that 
degenerates as t — > 00 or r — > t. 



References 

[1] Akhiezer, A., Akhiezer, I., Polovin, R., Sitenko, A., and Stepanov, K. Plasma 
electrodynamics. Vol. I: Linear theory, Vol. II: Non-linear theory and fluctuations. Pcrgamon 
Press, 1975 (Enlglish Edition). Translated by D. tcr Haar. 

[2] Caglioti, E., and Maffei, C. Time asymptotics for solutions of Vlasov-Poisson equation 
in a circle. J. Statist. Phys. 92, 1-2 (1998), 301-323. 



LANDAU DAMPING 



9 



[3] Degond, P. Spectral theory of the linearized Vlasov-Poisson equation. Trans. Amer. Math. 

Soc. 291 2 (1986), 435-453. 
[4] Glassey, R., and Schaeffer, J. On time decay rates in Landau damping. Comm. Partial 

Differential Equations 20, 3-4 (1995), 647-676. 
[5] Hwang, J.-H., and Velazquez, J. On the existence of exponentially decreasing solutions 

of the nonlinear landau damping problem. Preprint, 2008. 
[6] Landau, L. On the vibration of the electronic plasma. J. Phys. USSR 10 (1946), 25. English 

translation in JETP 16, 574. Reproduced in Collected papers of L.D. Landau, edited and with 

an introduction by D. tcr Haar, Pcrgamon Press, 1965, pp. 445-460; and in Men of Physics: 

L.D. Landau, Vol. 2, Pcrgamon Press, D. tcr Haar, ed. (1965). 
[7] Malmberg, J., Wharton, C, Gould, R., and O'Neil, T. Plasma wave echo experiment. 

Phys. Rev. Letters 20, 3 (1968), 95-97. 
[8] Maslov, V. P., and Fedoryuk, M. V. The linear theory of Landau damping. Mat. Sb. 

(N.S.) 127(169), 4 (1985), 445-475, 559. 
[9] Mouhot, C, and Villani, C. On the landau damping. Available online at 

http://arxiv.org/abs/0904.2760. Preprint, 2009. 
[10] Saenz, A. W. Long-time behavior of the electic potential and stability in the linearized 

Vlasov theory. J. Mathematical Phys. 6 (1965), 859-875. 
[11] Villani, C. Hypocoercivity. To appear in Mem. Amer. Math. Soc. 

[12] Villani, C. A review of mathematical topics in collisional kinetic theory. In Handbook of 
mathematical fluid dynamics, Vol. I. North-Holland, Amsterdam, 2002, pp. 71-305. 

[13] VLASOV, A. A. On the oscillation properties of an electron gas. Zh. Eksper. Teoret. Fiz. 8 
(1938), 291-318. 

Clement Mouhot 

University of Cambridge 
DAMTP, Centre for Mathematical Sciences 

WlLBERFORCE ROAD 

Cambridge CB3 0WA 
ENGLAND 
On leave from: 
ENS Paris & CNRS 
DMA, UMR CNRS 8553 

45 RUE d'Ulm 
F 75320 Paris cedex 05 
FRANCE 

E-MAIL: Clement.Mouhot@ens.fr 



Cedric Villani 



10 



C. MOUHOT AND C. VILLANI 



ENS Lyon & Institut Universitaire de France 
UMPA, UMR CNRS 5669 

46 ALLEE D'lTALIE 

69364 Lyon Cedex 07 
FRANCE 



E-MAIL: cvillani@umpa.ens-lyon.fr 



